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ABSTRACT 

We employ detailed photoionization models to infer the physical conditions of intrinsic narrow absorption 
line systems found in high resolution spectra of three quasars at z = 2.6-3.0. We focus on a family of intrinsic 
absorbers characterized by N V lines that are strong relative to the Lya lines. The inferred physical conditions 
are similar for the three intrinsic N V absorbers, with metallicities greater than 10 times the solar value (assum- 
ing a solar abundance pattern), and with high ionization parameters (logC/ ~ 0). Thus, we conclude that the 
unusual strength of the N V lines results from a combination of partial coverage, a high ionization state, and 
high metallicity. We consider whether dilution of the absorption lines by flux from the broad-emission line 
region can lead us to overestimate the metallicities and we find that this is an unlikely possibility. The high 
abundances that we infer are not surprising in the context of scenarios in which metal enrichment takes place 
very early on in massive galaxies. We estimate that the mass outflow rate in the absorbing gas (which is likely 
to have a filamentary structure) is less than a few M© yr _1 under the most optimistic assumptions, although it 
may be embedded in a much hotter, more massive outflow. 

Subject headings: quasars: absorption lines — quasars: individual (HE0130— 4021, Q1009+2956, 
HS 1700+6416) 



1. INTRODUCTION 

Intrinsic absorption lines seen in the spectra of quasars 
are thought to originate in ionized g as physically relate d 
with the quasar "central engine" (e.g.. lHamann et all 120041) . 
They include resonance lines that appear in the rest-frame 
UV band, such as N V, C IV and/or high ionization species 
that have absorption lines and edges in the soft X-ray band, 
such as O vii and O vm (e.g., lReynoldsl[l997t iGiustini efafl 
2010). The intrinsic UV absorption lines are usually 
blueshifted relative to the quasar emission lines and they 
are believed to arise from gas that is outflowing from the 
immediate vicinity of the central engine. Lines that are 
blueshifted by less than 5000 km s are conventionally re- 
garded as associated absorption lines (AALs). Many of 
these may i ndeed be connecte d to broad emission-line region 
(BELR) dGangulv et al.l200lh and to X-ray "warm" absorber s 
(iMathur. Elvis & Wilkes 111991 iBrandt. Laor & WiUslHOQOT) . 
These absorption lines thus trace the physical conditions and 
kinematics surrounding the central engine over a wide range 
of redshifts, and thus provide a diagnostic of quasar evolution. 

Narrow intrinsic absorption lines (NALs) typically have 
widths < 500 km s so that thei r doublets can be sepa - 
rated at high spectral resolution (see lHamann & Sabra 2004). 
This is in contrast to broad absorption lines (BALs) whose 
width can reach 30,000 km s" 1 dTrump et all 120061) . Intrin- 
sic NALs studies are complementary to studies of BALs, 
since they probe different regions in the quasar surround- 
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ings. Although the intrinsic NALs are more difficult to iden- 
tify, several practical considerations make the derivation of 
their physical conditions m ore straightforward than for BALs. 
(Hamann & Ferland 1999). Because NALs are often unsatu- 
rated and resolved, we can measure directly the NAL cover- 
age fractions and column densities of various ions. In BAL 
systems, doublet transitions are often self-blended, making 
such measurement more difficult. 

App roximately 60% of all quasars show evidence o f out- 
flows dMisawa et al.ll2007b IGanguly & Brothertonll2008l) . Al- 
though statistical methods can be used to determine this, par- 
tial coverage of doublets and multiplets and time variability 
analysis are the two most commonly used ways to determine 
decisively if a particular NAL system is intrinsic. These sig- 
natures are only seen in low ionization transitions of interven- 
ing absorbers, and only rarely in the case of small molecu- 
lar clouds that could be smaller than the projected size of the 
background quasar broad emission line region, (e.g., Jones 
et al. 2010, Ivanchik et al. 2010). Partial coverage is de- 
scribed by the coverage fraction Cf, which is the fraction of 
photons f rom the background source that p ass through the 
absorber (Barlow, Hamann, & Sargent 1997). It can be esti- 
mated using the residual flux ratio of resonanc e doublets (e.g., 
iBarlow & Sargenllll997t IGanguly et al.lll999l) . Variability of 
the absorption lines could be caused by transverse motion of 
the absorbing material or by chang es in its ionization state 
dHamannlll997l iMisawa et aD 120051) . Using UV spectra of 
z ^ 1.5 quasars observed at two different epoch s separated 
by 4- 10 years with the Hubble Space Telescope. IWise et al.1 
(2004) concluded that a minimum of 21% of the AALs are 
va riable. A similar conclu sion was reached for z ~ 2 quasars 
bv lNaravanan et al.1 (120041) . 

Early work investigating metal abundances using narrow 
associated absorptio n lines indicated s upersolar me tallici- 
ties, i.e.. Z>Zr., dPetitjean et al.1 Il994t iTripp. Lu & Savagd 
Il996t lHamannlll997l) . More recent work measures metallic- 
ities consistent with the previously determined values, and 
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TABLE 1 

Summary of Properties of Intrinsic NAL Systems a 
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Transitions Detected 


HE0130-4021 


3.030 


25.7 


48.36 


58.40 


2.9739 


4204 


3612-6084 


Lya, Ly/3, ClII,Nv, O VI 


Q 1009+2956 


2.644 


21.8 


48.49 


58.53 


2.6495 


-453 


3093^4623 


Lya, Ly/3, C III, S VI, N V, O VI 


HS 1700+6416 


2.722 


22.6 


48.98 


58.97 


2.7125 


764 


3726-6190 


Lya, Lyj3, C IV, N V, O VI 



" Adopted from Misawa et al. 1 2007). 

k The luminosity distance, assuming Q\ — 0.7, £2m — 0.3, h — 0.7. 

c The log of the bolometric luminosity (in erg s _1 ). See details in i|3.2| of the text. 

The log of the ionizing photon rate (in s~'). See details in i|3.2l of the text. 
e Absorption system redshift determined by finding the wavelength that divides the optical depth of the N V A 1239 line in half. 
^ Ejection velocity relative to emission lines. A positive sign denotes a blueshift. 
g Wavelength range covered by the Keck spectrum. 



in some cases, it ha s provided a more stringent constraint. 
iPapovich et alj (120001) analyzed az = 1 .207 AAL system to- 
wards Q0122+0388 with the help of photo ionization mod- 
els and found Z ~ 2Z Q . iFields et al. (2005) found that the 
metallicity must be at least 5Z Q in order to produce the 
high-ionization absorption lines (O VI, C IV, and N v) of the 
narrow -line Seyfert galaxy Mrk 1044. iGabel. Arav & Kiiril 

(2006) also found supersolar abundances by modelling the 
AALs of the quasar J2233-606, with [C/H], [O/H] « 0.5- 
0.9, and [N/H] ks 1.1 - 1.3. High metallicity up to 3Z W , is 
sugge sted in the study of quasar PG 0935+417 ( Hama nn et al 
2004). A study of the outflow from Mrk 279 by lArav et al 

(2007) has found an overabundance carbon, nitrogen, and 
oxygen relative to the Sun by factors of 2.2 ± 0.7, 3.5 ± 1.1, 
and 1.6 ±0.8, respectively. Moreover, photoionization mod- 
eling of associated NALs systems find metallicities ranging 
between a few and 20 times the Solar value (e.g., Petitjean et 
al. 1994, Tripp, Lu & Savage 1996, D'Odorico et al. 2004) . 
Such high metal abundances are also inf erred from ana l ysis of 
the broad emission lines. For example. lHamann et alJ (2002) 
find that a-elements, such as C and O, are overabundant in 
the broad-emission line regions of quasars by a factor of ~ 3 
relative to the Sun. 

Previous studies indicate that the AALs and BALs are lo- 
cated within the scale of the host galaxy and suggest that the 
absorbers are physically associated with and may originate 
from the quasar outflow. For example, the strengths of the 
excited-state Si II* AALs in the high resolution spectrum of 
3C 191 indicate a distance of 28 kpc from the q uasar assuming 
the gas is photonionized (Ha mann et al.1 120011) . By studying 
the me tastable excited levels of Si II and Si II*, Du nn et alJ 
(2010) determines the distance to the outflow for attenuated 
and unattenuated SEDs of ~ 6 kpc or ~ 17 kpc respectively. 
Variability analysis of a C IV NAL complex towards quasar 
HS 1603+3820 constrains an upper l imit on the distance of 
r < 6 kpc from the continuum source (Misa wa et al.l l2005). 

Although it is generally agreed upon that intrinsic absorp- 
tion lines are formed in an outflow from the central engine, 
there is no consensus on the geo metry or the d r iving force of 
such an outflow. In the model of Murray et al.l (1 19951). w hich 
was developed further by Prog a. Stone & Kallman (200(1), the 
wind is launched from the entire surface of the accretion disk, 
with many denser, absorbing filaments distributed throughout. 
In this picture, the BALs appear at low latitudes (small angles 
from the plane of the disk) when the line of sight passes di- 
rectly through the dense, high velocity wind. The NALs could 
form at higher latitudes when the line of s ight intercep ts a 
dense filament in the hot, low-density flow. Elvis (2000|), in- 



stead, proposed a funnel-shaped thin shell outflow model, in 
which the wind rises vertically from a narrow range of radii 
on a disk in the BELR. The high-ionization UV NALs and 
the X-ray "warm absorbers" are seen when the disk is viewed 
at a low latitude, through the wind. When viewed along the 
fast outflow, the full range of velocity is seen in absorption 
with a large total column density, giving rise to a BAL. Re- 
cent work indicates that the column densities measured in 
the X-ray spectra of quasars ho sting intrinsi c NALs are con- 
siderably lower than what the lElvisI (120001) scenario posits 



2009), which suggests a 



(Misa wa et all 12008; Chartas et al.1 . 

preference on the model of lMurrav et al.l d 19951) . There is also 
a suggestion that in some objects, the outflows are directed in 
the polar direction, parallel to the axis of the accretion disk 
(see, for example. iGhosh & Punslvl l2007t IZhou et al J 120061 
iBrotherton et al.ll2006l) . In the context of such a picture, and 
by analogy with equatorial wind models, the absorbers giving 
rise to NALs could be detached from the dense part of the out- 
flow and become visible along different lines of sight. Previ- 
ous studies have shown that the absorbing clumps originate at 
similar physical locations and are driven by radiative acceler- 
ation. This is supported by the explanation of a double trough 
in the C IV BAL by the "line-locking" effect (e.g., Ar a v et al . 
1995). A striking case can be seen in Srian and et al.l (2002), 
which reports an observation of a highly structured flow in 
w hich distinct compon ents have similar velocity separations. 

Misa wa et al.l d2007l) identified two families of intrinsic ab- 
sorption lines among 39 intrinsic NAL systems. The "strong 
N V" family is characterized by strong, partially covered N V, 
relatively weak Lya (less than twice the equivalent width of 
N v), occasionally detected C IV, and O VI lines (sometimes 
even stronger than the N V lines). The "strong C IV" family, 
on the other hand, is characterized by strong, partially cov- 
ered C IV doublets, strong, usually black Lya lines, and rela- 
tively weak or undetected N V. A possible thir d family, char- 
acterized by strong O VI lines is discussed by Ganguly et al. 
(12010. in pr eparation). They primarily used partial covering 
to diagnose the intrinsic nature of the systems, i.e., they found 
that the members of the doublet could not be fit simultane- 
ously assuming full coverage. Depending on the geometries 
and locations of the continuum source and broad emission 
line regions (and empirically on their positions relative to the 
emission lines) the ab sorber can partially cover either or both 
(Ganguly et al. 1999). There may be overlap between the C IV 
and O VI categories as some absorption systems can possess 
both strong C IV and O VI lines. Further observations cov- 
ering both of these lines are needed to investigate this issue. 
These different types of NALs are of interest because they 
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may allow us to probe different regions of the outflow. These 
three families of NALs may represent different lines of sight 
through the outflow, which is also suggested by the relations 
between the properties of UV NALs and the X-ray prope rties 
of the quasars that display them (e.g.,[Chartas et al. 20091). 

The question we address in this paper is the origin of the 
intrinsic N V NALs. To this end we construct photoionization 
models for the strong N V absorption systems in the spectra of 
three r adio-quiet qua s ars fr om the HIRES/Keck sample stud- 
ied by iMisawa et alJ (120071) . This sample contains 37 opti- 
cally bright quasars at z = 2 — 4. We choose these particu- 
lar three quasars because their spectra exhibit characteristics 
of the "strong N v" family. In addition, these three systems 
offer many observational constraints because many ions are 
covered in their spectra, and the Nv absorption lines have 
relatively simple profiles. Similar systems, containing multi- 
ple absorption components spread over thousands of km s _1 
have been reported in the spec tra of RX J123 0.8+0115 
(iGaneulv et alJl2003l) and 3C 351 (lYuan et al.H2002h . 

In Section|2] we describe the absorption profiles of the three 
intrinsic Nv NAL systems. In Section [3] we introduce our 
method for modell ing the three sy stems using the Cloudy pho- 
toionization code (Ferland 2006). The modelling results, in 
the form of constraints on metallicity, ionization parameter, 
and volume number density are presented in Section [4] In 
Section[6] we discuss possible interpretations of our results on 
the location of the absorbers in the quasar winds. We present 
a summary and conclusion in Section [7] The cosmology we 
use in this paper is £l\ = 0.7, £2m = 0.3, h = 0.7, which leads 
to the luminosity distances listed in TableQ] 

2. DATA 

The Keck/HI RES spectra of our th ree intrinsic N V systems 
are described in Mis awa et alJ (120071) . The spectral resolution 
is R = 37,500, or ~ 7 km s _1 . Table Q] summarizes the basic 
properties of the three absorption systems and of the quasars 
that host them. The transitions listed are detected at a 5(7 con- 
fidence level at the NAL redshift. The sign of the velocity of 
a system is taken to be positive, if the line is blueshifted, i.e., 
if the gas appears to be outflowing relative to t he quasar; this 
is the opposite convention from that adopted in Misa wa~et alJ 
(120071) . In the s ame table, we also l ist the 4400 A flux densi- 
ties (taken from|Misawa et al. 2007), as well as the bolometric 
luminosities and ionizing photon rates (obtained as described 
in ®. 

Figures[T][2] and|3]present, for each system, absorption pro- 
files for transitions which are used as modelling constraints. 
An example of the best models we have found (discussed in 
later sections) is also shown in each figure. The velocity of an 
entire system is defined by the optical depth-weighted center 
of the strongest member of the N V doublet. We describe each 
system below. 

: HE0130-4021 (z abs = 2.973915; Fig. [Qi: The system is 
found ~ 4000 km s blueward of the quasar emission 
redshift, just within the associated region. This sys- 
tem is kinematically more complex than the other two, 
with 9 components in each N V line, spreading over 
~ 500 km s , All but one of these components have 
small coverage fractions (Cf < 0.5; see Table |2). The 
strong features in the Lya panel of FigureQ]at ~ —280, 
50, 110, and 210 km s cannot be Lya components 
because the corresponding Lyj3 components are not de- 
tected. Thus only portions of the Lya profile can be 



used as constraints, but it is clear that the absorption in 
Lya is quite weak. The O VI doublet appears to be de- 
tected in many of the same components as N V, though 
inconsistencies between Ovi A 1032 and Ovi A 1038 
betray a number of blends. For this system the align- 
ment of kinematic components of different transitions 
is convincing and it seems at least some of the Ovi ab- 
sorption is real. We can use the observed absorption at 
the positions of O VI AA 1032, 1038 as an upper limit. 

: Q1009+2956 (z abs = 2.649510; Fig. |2): This system is 
~ 450 km s redward of the quasar emission red- 
shift. The N V AA1239, 1242 doublet of this system 
exhibits partial coverage in each of the three blended 
kinematic (Voigt) components required to fit its asym- 
metric profile, with the strongest component having 
C f = 0.70 ± 0.06 (Fig. and Table 0. Lya and Lyj3 
are detected, but are unusually weak compared to N V. 
The ratio of Lya to Lyj3 is reasonable, but differences 
in the shapes of these lines, particularly the presence 
of a wing on the red side of the Lyj3 profile, suggest a 
possible blend with another line. We will thus weight 
the Lya line more heavily as a constraint. Narrow 
Cm A 977 and S VI A933 lines appear to be detected, 
albeit with some uncertainty because of their location 
in the Lya forest. The red member of S VI AA933,945 
doublet is heavily blended and thus cannot be used to 
judge whether this doublet is truly detected or not. The 
feature at the position of Si III A 1207 is not aligned 
with the C III A977 and N V A 1239 profiles and is too 
broad, so we believe that we only have an upper limit 
on Si III A 1207 absorption from this system. There 
is absorption at the positions of both members of the 
Ovi AA 1032, 1038 doublet, however the minima are 
not aligned with each other or with the N V profile. The 
Ovi A 1032 is too broad relative to Ovi A 1038, and 
both transitions are clearly affected by blends. The lim- 
its on O VI for this system are not restrictive constraints 
for our models. 

: HS1700+6416 (z abs = 2.7125; Fig. [3): This system is ~ 
750 km s _1 blueward of the quasar emission redshift. 
The Nv AA 1239, 1242 and Civ AA 1548, 1550 dou- 
blets are detected in this system, as well as Lya and 
Ly/3 (Fig. [3]>. The same s ystem was presented in 
Barlow, Hamann, & Sar gentl (119971) . who found it to 
vary significantly on a time scale of ^6.5 months in 
its rest-frame. The central absorption is detected at 
z~ 2.7125 (u re i ~0kms _l ) and a weaker component is 
found at ~ 325 km s . The central component shows 
an asymmetry, indicating the need for two Voigt com- 
ponents for an adequate fit. The feature at ~ 70 km s _1 
in the Lya panel cannot be Lya belonging to this ab- 
sorption system, since no counterpart is detected in the 
Lyj3 panel. All three components of the Nv profile 
show partial coverage, with coverage fractions between 
0.33 and 0.43 (Table 01. The Lyj3 component at km 
s is too strong relative to Lya, indicating a blend that 
affects the Lyj3 so that it can only be used as an upper 
limit. If saturated (and affected by partial coverage), 
Ly/3 could have an equivalent width equal to that of 
Lya, but it can never have a larger equivalent width. 
The O VI doublet lies in the Lya forest, such that it is 
hard to evaluate whether it is detected for this system or 
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TABLE 2 

Results of Decomposition of N v Doublet Profiles a 



Kinematic u e j b u re i c b c 

Quasar Component (kms~') (kms~') logfA'Nv/cm -2 ) d (kms -1 ) Cf f 



HE0130-4021 



Q1009+2956 



HS 1700+6416 



4324 
4275 
4236 
4209 
4084 
4020 
3891 
3952 
4210 

-453 
-455 
-461 

771 

766 
450 



-120 
-71 
-32 
-5 
+ 121 
+ 184 
+313 
+252 
-6 

-5 
+2 
+7 

-9 
+2 
+314 



14.18 + 0.04 
14.2 + 0.1 
13.9 + 0.1 

14.8 + 0.3 

15.9 + 0.5 
15.0 + 0.3 
15.2 + 0.6 

14.62 + 0.06 
13.51+0.01 

13.70 + 0.29 
14.79+1.48 

14.63 + 0.28 

14.18 + 0.05 
14.62 + 0.04 
13.8 + 0.1 



16.6 + 0.9 
20 + 8 
13 + 2 
32 + 5 
28 + 5 
17 + 2 
6+1 
34 + 3 

12.6 + 0.3 

14.3 + 1.7 
3.40+1.90 

13.0 + 2.1 

11.0 + 0.5 
25.6 + 0.5 

18.4 + 0.8 



0.43 + 0.03 
0.16 + 0.04 
0.49 + 0.08 
0.13 + 0.03 
0.10 + 0.02 
0.13 + 0.03 
0.37 + 0.02 
0.21+0.02 
1.00 g 

0.62 + 0.11 
0.70 + 0.06 
0.22 + 0.05 

0.43 + 0.02 
0.33 + 0.02 
36+ 07 



a The errors are a combination of the continuum fitting errors and Voigt profile fitting errors, 
k Ejection velocity relative to the emission lines. A positive sign denotes a blueshift. 

c Velocity offset relative to the redshift of the entire absorption line complex. The sign of the velocity is set to be positive, 
if the lines are redshifted from the system center. 

The N V column density inferred from the Voigt profile fit. 
c Doppler parameter. 
^ Effective coverage fraction. 

g There is no error bar for this entry because our fitting yields Cj > 1, which is unphysical. Thus, we set Cr — 1.00 when 
we fit this component. 



not. If all of the detected absorption, either at km s 1 
or at 325 km s , is due to Ovi, the absorber would 
have a much larger coverage fraction than the N V, C IV, 
or Lya absorbers at these same velocities. In our mod- 
elling we will take the absorption lines coincident with 
O VI as an upper limit, but we note that the minima of 
the absorption at ~ 325 km s _1 in both members of the 
Ovi doublet are at the same velocity as those in the 
N V. Because of this, we will later consider whether vi- 
able single phase models can also explain this possible 
O VI absorption. 

Transitions appearing in Figures Q] [2] and [3] but not de- 
scribed above, provide limits for model constraints, either be- 
cause they are not detected or because they thought to be af- 
fected by blends. 

3. PHOTOIONIZATION MODELING 

Our calculations were performed with version 07.02.01 of 
Cloudy, last described by iFerlaridl d!998f) . Our modeling as- 
sumes a series of plane-parallel slabs of gas (clouds) exposed 
to the ionizing continuum from the central engine. The gas 
within each absorbing cloud is assumed to have a uniform 
density, metallicity and abundance pattern. In most c ases, we 
assume a solar abundance pattern (Holweger 2001), so that 
for every element i, A^/A;.© = Z/Zq where A,- ~ tii/ «h- When 
an adequate fit cannot be achieved using the solar abundance 
pattern, we consider deviations. The intensity of the ionizing 
continuum is parameterized by the ionization parameter, 



U 



«H 



1 



\%r L 



*2L P 



cn H J £l £ 



■de, 



(1) 



where r is the distance of the illuminated face of the cloud 
from the continuum source, £\ is the lowest energy required to 



photoionize the gas, i.e., £i = 1 Ry, and £2 is the high energy 
cutoff of the SED. In our model, we take e 2 = 7.354 x 10 6 Ry. 
Each of our models are specified by U and tin, which (given 
the quasar luminosity and spectral shape) corresponds to a 
given r. We also assume that all the ions within a cloud have 
the same coverage fraction, and that the gas cloud is in a state 
of thermal equilibrium, so it is parameterized by a single elec- 
tron temperature, T e . The assumption of the same coverage 
fraction for different lines is likely to be valid for those from 
ions with similar ionization states, particularly C IV and N V 
which are the most important constraints for out models. Lya 
absorption may also arise from additional regions, but in any 
case will provide a lower limit on metallicity. In our favored 
models, the clouds are optically thin to the incident continuum 
so that the incident SED does not change after the ionizing 
photons pass through an absorbing cloud. 

We have measured the Doppler b parameters, column densi- 
ties and coverage fractions of the N V doublets, as we describe 
below, and optimized on these values (i.e., we required that 
Cloudy models produce them). We choose Nv as the transi- 
tion on which we optimize because it is located in a relatively 
"clean" spectral region where there are few blends. This al- 
lows us to determine its column density, Doppler parameter, 
and co verage fraction via Voigt profile fit ting. The code AU- 
TOVP (iDave. Hernquist & Weinberglll997l) is used to derive 
an initial solution, and then MINFIT is used to determine the 
minimum number of components that produce an adequate 
fit. The goal of the modeling exercise is to reproduce the ob- 
served absorption profiles for all other ions by adjusting three 
physical parameters: metallicity, Z/Z Q , ionization parameter, 
U, and hydrogen number density, «h- 

3.1. Voigt Profile Fitting of the Nw Doublets 

The approach of the MINFIT code is to first "overfit" the 
system using many Voigt components and then to reject com- 
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FIG. 1. — Plots of absorption lines from all the transitions used as con- 
straints for the v e j = 4204 km s (z a (, s = 3.973915) absorption system in the 
spectrum of HE0130— 4021. Transitions are ordered by increasing ionization 
potential, except that Lya and Lyj3 (Hi 1026) are always the first. The zero- 
point of the velocity scale is set by the bisector of the strongest blend of the 
Nv line (the bisector divides the optical depth of the blend in half). Over- 
plotted (red) in each panel is the best theoretical profile we obtained by ad- 
justing the photoionization model parameters (see Section|4]for details). We 
mark each kinematic (Voigt) component using a tick mark; the component 
numbers are printed in the Nv A 1239 panel. The corrsponding componets 
for different transitions have the same partial coverage factor (see SectionfJJ- 
The corresponding components for different transitions have the same partial 
coverage factor (see SectionfJJ, which are listed in TablefS] 

ponents that do not improve the fits at a confidence level above 
95%. This fitting technique has been used extensively in stud- 
ies of intervening Mg II systems (e.g., Ding et al. 2003; Ding, 
Charlton & Churchill 2005; Lynch, Charlton & Kim 2006). 
Table|2]lists the resulting fitting parameters. It is essential to 
include the coverage fractions as free parameters in the fit- 
ting process for these systems since Cf = 1 does not provide 
acceptable fits to these doublets within the observed uncer- 
tainties (Fig.gl Misawa et al. 2007). 

3.2. Specification of the Ionizing Continuum 

Although multi-band photometry and/or spectra of these 
three quasars are not available over the entire spectral range, 
some observational data are available to place constraints on 
the input SED. We consider three possible quasar SEDs. 



1. A typical quasar SED from lElvis eTail (119941) . This is 
a mean SED determined from 29 moderately luminous 
quasars. 



2. A bro ken power-law SED from Mat hews & Ferlandl 
(119871) . This is a mean continuum determined by com- 
bining direct observations, without distinguishing be- 
tween radio quiet and radio loud quasars. 



FIG. 2.— Same as Fig.[T]but for the u eJ = -453 km s -1 (z abs = 2.649510) 
absorption system in the spectrum of Q1009+2956. 
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FIG. 3.— Same as Fig.Q]but for the u cj = 764 kms ' (z abs = 2.7125) 
absorption system in the spectrum of HS1700+6416. 
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FIG. 4.— Voigt profile fits to the N V AA1239, 1242 doublet in the u cj = 
764 km s -1 (z abs = 2.7125) system in the spectrum of HS1700+6416. The 
velocity zero point is defined as in Fig. [T] The fitting process minimizes the 
number of kinematic components needed for an adequate representation of 
the profile. Vertical ticks mark the component centers. In the upper pair of 
panels the coverage factors of all components are fixed to be unity; in the 
lower pair of panels, the coverage factors of all components are treated as 
free parameters. Black dots with error bars represent the value of (1 — Cf) 
for each component. Clearly, Cf = 1 does not provide an adequate fit, so that 
partial coverage is required. 

3. A multi-component continuum with either a small, a 
medium, or a big blue bump (BB) (see the Cloudy 
document, Hazy § 6.2, for details), parameterized by 
fc!T RR = lQeV, 100 eV and 30 eV, respectively 
(Casebeer, Leighly, & Baron 2006). By adjusting the 
parameters of this prescription we produce three dif- 
ferent strengths of the BB, which we list in Table [3] 
("small","medium", and "big"). For each BB strength 
we examine three different values of Oo X (denoted as I, 
II, and III in TableO, which gives us nine different SED 
shapes for this model. 

We utilize three parameters to constrain the SED: a a , Ox, 
and Oqx, which are the UV/optical spectral slope, X-ray 
band spectral slope and optical-X-ray spectral slope, respec- 
tively (assuming that f v v +a ). The last is defined as 
a ox = 0.3838 log [£ v (2 keV) /£ v (2500 A)] , where £ v {2 keV) 
and £ v (2500 A) are the monochromatic luminosities at 2 keV 
and at 2500 A dTananbaum et al.ll 19791) . 

The low resolution (R ~ 2000) UV/optical spectrum of 
HS 1700+6416 is available from the D ata Release 5 (DR5) 
quasar catalog ( Schneid er et all 120071) of the Sloan Digital 
Sky Survey (SDSS; York et al. 2000). We fitted the underly- 



TABLE 3 

Spectral Indices of Different SED Models 
Compared to Observed Values 



SED a 


fv 






Observed, HS 1700+64 16 


-0.47 


-1.00 


-1.87 


Elvis 1994 


-0.91 


-0.93 


-1.37 


Mathews 1987 


-0.50 


-0.70 


-1.41 


Small BB 1 


-1.24 


-0.99 


-1.40 


Small BB II 


-1.27 


-0.99 


-1.60 


Small BB III 


-1.27 


-0.99 


-1.90 


Medium BB I 


-0.54 


-1.00 


-1.40 


Medium BB II 


-0.56 


-1.00 


-1.60 


Medium BB III 


-0.56 


-1.00 


-1.90 


Big BB I 


-0.49 


-1.82 


-1.40 


Big BB II 


-0.50 


-2.30 


-1.60 


Big BB III 


-0.51 


-3.06 


-1.90 



' The Roman numerals 1,11,111 represent different values of a m 




16 IS 
log[v (Hz)] 

FIG. 5. — Shapes of two SEDs we used as input to our photoionization 
models. The solid line shows the SED used for HS 1700+6416 while the 
dotted line shows the SED used for the other two quasars. These SEDs are 
normalized to an ionization parameter \ogU = 0.1 for a hydrogen density of 



log(fi H /cm 



Thus, they have the same photon luminosity between 1 



and 7.354 x 10 Ry. Also labelled are positions of important transitions - 
Lyce, Nv A A 1239, 1242 and Civ XX 1548, 1550. We also mark the photon 
frequencies needed to ionize N+ 3 to N+ 4 , C+ 2 to C+ 3 , and H° to H+ . We also 
mark the frequencies corresponding to 2500 A and 2 keV since flux densities 
at these two frequencies are used to define a ox . 

ing continuum using a power-law, with a s mall Balmer bump 
and the iron emission forest superimposed (IVanden Berk et all 
2010, in preparation), and obtained a Q = —0.47 (see TableO. 
This object was also observed by Chandra and XMM-Newton. 
The X-ray spectrum was adequately fit with a power-law in- 
dex ~ — 1 (oc x = —1.2 ±0.2 for the Chandra spectrum and 
a x = -1.1 ±0.2 for the XMM-Newton spectrum; Misawa et 
al. 2008). Thus we adopted a x = - 1 . 

To find a ox we used the empirical correlation between it and 
the 2500 A quasar luminosity ( Vignali 2003^ IStrateva et alj 
[2005t ISteffenetaL|[2q06b. W e adopt the a ox -^ v (2500 A) re- 
lation by lSteffen et all (120061) : 

a ox = ( -0.1 37 ± 0.008) log^ v (2500 A) + (2.638 ±0.240). 

The 2500 A luminosity is calculated as ^ v (2500 A) = 
4-nd^fv (2500 A), in which d^ is the luminosity distance and 
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/ v (2500 A) is the rest-frame flux, derived from an extrap- 
olation of the flux den sity at 4400 A dMisawa et alj 12007b 
assuming f v oc V - QM dVanden Berk et alj l200ll) . Thus, for 
HS 1700+64 16 we find = -1.87 

The three spectral indices for HS 1700+6416, from spectro- 
scopic observations, are compared to those for possible quasar 
SEDs (Table [3). We find that the Medium BB III SED best 
matches the observational constraints thus, we use this SED 
as the input ionizing flux for HS 1700+6416. This SED is 
plotted as a solid line in Fig. [5] 

We measured a of the other two quasars from the low- 
resolution spectra o f Q1009+2956 dBurles & Tvtleri [19981) 
and HE0 130-4021 dOsmer. Porter & Greenl 11994 . respec- 
tively. These UV/optical power-law indices are determined 
by connecting two empirical "line-free" regions redward of 
the Lya emission line and are ~ —0.25 for both quasars. The 
2500 A luminosities for these two quasars are determined in 
the same manner as for HS 1700+6416 and their a ox values 
are both ~ —1.8. There have not been any X-ray spectral 
observations of Q1009+2956 or HE0130-4021 so, we as- 
sume they have the same a x as HS 1700+6416. Generally, 
the photoionization cross section decreases rapidly with pho- 
ton energy so the ionization state is not sensitive to the X-ray 
spectral slope. We then also use the Medium BB III multi- 
component SED for both of these quasars, but we adjust the 
optical slope to —0.25 to match the observations. This SED 
is shown as a dotted line in Fig. [5] 

We have integrated the above SEDs to compute the bolo- 
metric luminosities and ionizing (E > 1 Ry) photon rates for 
the three quasars (Lbol an d respectively). We first normal- 
ized them to match the 4400 A luminosity density of each 
quasar (using the flux density and luminosity distance from 
Table [T). The resulting values are listed in Table [T] they will 
be useful for our estimates of the properties of the absorbing 
gas in Sj6]2. 



3.3. Comparison of Cloudy Model Results to the Data and 
the Role of the Coverage Fraction 

Cloudy calculates the electron temperature based upon ther- 
mal and ionization balance. Using this temperature, we cal- 
culate the microturbulence velocity for N V based on its ob- 
served total Doppler b parameter assuming that both of them 
follow Gaussian distribution. We then calculate the b pa- 
rameters for all other elements by combining this turbulent b 
with the thermal b, calculated from the atomic weight and the 
electron temperature. Using these b parameters and the col- 
umn densities of each ion, output by Cloudy, we synthesize 
noiseless model spectra, convolving with a Gaussian function 
which represents the instrumental profile. We compare the 
theoretical profiles with the observed ones using a % 2 test in 
combination with visual inspection. The wavelength range 
within which x 2 is evaluated is carefully chosen because even 
a single pixel affected by a blend or an instrumental a artifact 
can dominate the % 2 values. 

We apply a grid method in (Z/Z & )-U-nn space to search 
for acceptable solutions by comparing model profiles to the 
observed absorption lines. In order to better visualize the 
models and compare them to the observed spectra, we take 
constant-«H slices and examine combinations of (Z/Z Q ) and 
U in this 2-dimensional parameter space. The ranges of 
(Z/Zq), U, and «h we explore vary depending on the spe- 



cific system, but typically they span 1 < log (Z/Zq) < 30, 
-1.5 < logt/ < 1.5, and 2 < log(n H /cm- 3 ) < 14. The ini- 
tial increments are Alog (Z/Zq) = 0.1, Alogt/ = 0.1, and 
Alog(«H/cm~ 3 ) = 2. Finer grids are applied only if neces- 
sary. We assume that different kinematic (Voigt) components 
represent different parcels of gas, which we model separately. 
However, if kinematic components are blended together, we 
start by assuming that they are all described by the same 
model parameters. We refer to such models as the "initial" 
or "preliminary" models. After finding the best preliminary 
solution, we vary (Z/Zq) and U for each kinematic compo- 
nent within 0.2 dex of their preliminary values to seek better 
solutions with smaller % 2 values. 

In the process of producing synthetic line profiles for com- 
parison to the observed spectra, we make use of the coverage 
fraction inferred from the relative strengths of the lines in the 
N V doublet (see Table [2]). Since Cf represents the fraction of 
photons from the background source that pass through the ab- 
sorber, we effectively dilute all the absorption lines of a given 
system by a factor equal to Cf. We have to assume that the 
same value of Cf applies to all transitions for lack of additional 
information. This n eed not be the case in general, however; 
Misa wa et al.l ((2007) do find cases where different resonance 
doublets in the same system yield different values of Cf, for 
example. In one of our three quasars, HS 1700+6416, we do 
have separate measurements of Cf for the N V and C IV dou- 
blets but they agree with each other within uncertainties. We 
return to the issue of the covering factor in Section where 
we discuss the possibility of the absorbers covering different 
fractions of the continuum and broad-emission line sources. 

4. PHOTOIONIZATION MODEL RESULTS 

4.1. HE01 30-4021 

This absorption complex can be described by nine kine- 
matic components (Fig. Q] and Table |2). These components 
are blended so we model them by first assuming that they 
are described by the same model parameters. Si IV XI 394 
is not detected, which yields a lower limit on logf/. Al- 
though Ovi XX 1032, 1038 and Cm X911 do not provide 
strong constraints, we require that they are not over-produced 
by our models. The Lya line provides a strong constraint for 
this system. Because the saturated absorption troughs around 
50 km s and 100 km s _1 are not seen in Lyj3, they can- 
not be due to Lyman absorption from this system. These two 
absorption features actually divide the whole Lya absorption 
system into a red part (positive velocity) and a blue part (neg- 
ative velocity), which we can model separately. The blue part 
comprises Components 1-4 and 9 while the red part consists 
of Components 5-8. Although this system comprises many 
kinematics components, we can still obtain robust results for 
the coverage fraction of each component. This is because we 
fit the blue and the red members doublet simultaneously so 
that only features appearing in both of doublets are contribute 
to the determination of Cf. The resulting relative uncertainties 
(listed in Table [2]) are typically of order 20%, which includes 
errors in profile fitting as well as errors in continuum fitting. 

We start by seeking models that can reproduce the observed 
Lya absorption complex. For either the red or the blue part 
of the profile, we can find a series of solutions for log (Z/Zq) 

5 In our initial simulations we tried to reproduce the absorption line 
strengths with log (Z/Zq) < 1. However, we found that acceptable models 
require relatively high metalicities, which led us to adopt this range of values 
in our grid search. 
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FIG. 6. — The log (Z/Zg/j-logJ/ parameter slice at log fin = 8, showing the tracks of acceptable preliminary models for the absorption system in the spectrum 
of HE0130— 4021. The left panel shows the tracks for the blue group (kinematic components 1—4, and 9). The right panel shows the tracks for the red group 
(kinematic components 5-8). Black dots are the cases we covered in our modelling. The dashed line represents the best-fit track for Lya while the dotted line 
represents the best-fit track for O VI. Arrows indicate the directions in which a specific theoretical absorption profile becomes stronger, so on the "tail" side, the 
absorption line is under-produced while on the "head" side, it is over-produced. 



and log U on a constant-/^ slice. Fig.[6]illustrates the method- 
ology we apply to all systems: it shows the acceptable models 
for Lya for log(«H/cm~ 3 ) = 8 and a range of log (Z/Z Q ) and 
logi/. The highest and lowest possible values of log (Z/Z Q ) 
and logf/ for different values of log(«H/cm -3 ) are tabulated 
in Table [4] The blue and red parts of the profile yield dif- 
ferent sets of acceptable model parameters. Fig. [6] presents 
constraints on models based on the observed Ovi A 1032 
and Si IV A 1394 lines. The effect of changing logf/ and 
log (Z/Z©) on the model profiles of Si IV A 1394 and Lya 
is illustrated in Figures [7a] and I7E1 The Ovi A 1032 pro- 
file is very insensitive to metallicity as well as to log U when 
logf/ > —1.5. Below this value, the Ovi A 1032 absorption 
decreases gradually. Because the O VI doublet is in the Lya 
forest, we use it as an upper limit on the O VI absorption that 
a model should produce. However, at logf/ ~ —1.5, the red 
parts of O VI 1032,1038 can be adequately fit. Therefore, the 
intersection between the acceptable models for the red part 
of Lya and the acceptable models for Ovi 1032 happens at 
logf/ ~ —1.5 and log (Z/Z Q ) ~ 1.7. This model also fits the 
observed C III A977 absorption and it does not overproduce 
Si IV A 1394. For the blue part of the system, within the range 
of models that are consistent with the Lya and SilV 1394, the 
O VI A 1032 is never overproduced. Since the O VI is taken as 
an upper limit, these models are all considered acceptable, as 
listed in Table @] 

The synthesized line profiles are consistent with the data 
for values of log(nn/ cm 3 ) between 2 and 14. The ranges of 



log (Z/Zq) and log U are slightly different for different log «h 
values (Table |U). An example of an acceptable model is pre- 
sented in Fig.Q] 

There are still discrepancies between the model and the ob- 
served spectrum of Lya. Perhaps the Lya line is blended 
with another transition because the corresponding absorp- 
tion trough is not seen in Lyj3. The v ~ 300kms _1 ab- 
sorption trough in the Lya window is also underproduced, 
which seems to suggest an additional weak component with 
a small Doppler parameter. Most of the absorption in O VI 
is underproduced by this model, with the exception of that at 
V ~ km s . The absorption at this velocity is mainly con- 
tributed by Component 9, which is the only component with 
a coverage fraction of unity. We have assumed that the O VI 
absorption is only an upper limit, but if it is real, a separate 
phase would be required that has a larger coverage fraction. 
In the other two systems O VI is si milarly underproduced; we 
discuss this further in Section I4~4l 

We have attempted to optimise the model parameters for 
each kinematic component separately but we found that this 
exercise did not yield useful results. We found a large num- 
ber of "pockets" in paramneter space, spanning a wide range 
of parameter values, where an acceptable fit can be achieved 
but there is no strong reason to prefer one particular "pocket" 
over the others. We attribute this outcome to the fact that the 
Lyman series lines that we use as constraints are blended with 
other lines and the profile of C III has a relatively low signal- 
to-noise ratio. Therefore, Table |4]provides the final summary 
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TABLE 4 

Summary of Preliminary Fitting Results 



Jvinematic 








Abundance 


Quasar Component 


log (Z/Z ) 


logf/ 


log(«H/cm 3 ) 


Pattern 


HE0130-4021 b 1-4,9 


1.6 


-0.6 to 


2 


solar 




1.6-1.7 


-1.4 to 


8 


solar 




1.5-1.6 


-1.4 to -0.3 


10 


solar 




1.5 


-1.3 


12 


solar 




1.3 


-1.0 


14 


solar 


5-8 


> 1.7 


> -0.6 


2 


solar 




1.9 


-0.9 to -0.2 


8 


solar 




> 1.8 


-0.8 to 0.4 


10 


solar 




1.6 


-0.8 


12 


solar 




1.8 


-0.6 


14 


solar 


Q1009+2956 1-3 


2.5 


0.2 


2-8 


solar, [S/H]< 1.5 




2.5 


0.3 


9 


solar, [S/H]< 1.5 




2.4 


0.2 


10 


solar, [S/H]< 1.4 




2.3 


0.3 


11-12 


solar, [S/H]< 1.5 




2.3 


0.7 


12 


solar, [C/H]~ 2.8 




2.5 


0.4 


13 


solar, [S/H]< 1.5 


HS1700+6416 a 1-2 


1.4 


0.5 


2-8 


solar 




1.3 


0.4 


9 


solar 




1.3 


0.5 


10 


solar 


3 


1.7 


1.0 


2-8 


solar 



NOTE. — Here we assume that blended kinematic components, as indicated in the second column of the table, are 
described by the same combination of model parameters. 

d The results presented in this table are also a summary of the final results for this quasar. See the discussion in jj4l l 
of the text. 



TABLE 5 

Best model parameters for the Q1009+2956 absorption 

SYSTEMS 



Kinem. 


log 




log 


Abundance 


r a 


Comp. 


(Z/Z ) 


log U 


(« H /cm~ 3 ) 


Pattern 


(PC) 


1 


2.4 


0.0 


2-8 


solar, [S/H]=1.3 


10- 


2 


2.6 


0.3 




solar, [S/H]=1.5 


10,000 


3 


2.3 


0.0 




solar, [S/H]=1.4 




1 


2.3 


0.1 


9 


solar, [S/H]=1.3 


-0.2 


2 


2.5 


0.3 




solar, [S/H]=1.5 




3 


2.7 


0.1 




solar, [S/H]=1.7 




1 


2.4 


0.0 


10 


solar, [S/H]=1.4 


-0.8 


2 


2.5 


0.3 




solar, [S/H]=1.5 




3 


2.3 


0.0 




solar, [S/H]=1.3 




1 


2.2 


0.2 


11-12 


solar, [S/H]=1.2 


0.04- 


2 


2.4 


0.3 




solar, [S/H]=1.4 


0.2 


3 


2.4 


0.5 




solar, [S/H]=1.4 




1 


2.2 


0.5 


12 


solar, [C/H]=2.7 


~ 0.04 


2 


2.4 


0.8 




solar, [C/H]=2.9 




3 


2.5 


0.5 




solar, [C/H]=3.0 




1 


2.7 


0.3 


13 


solar, [S/H]=1.7 


~ 0.004 


2 


2.6 


0.5 




solar, [S/H]=1.6 




3 


2.6 


0.4 




solar, [S/H]=1.6 





The distance of the absorbing gas from the ionizing source for this combination of U 
and «h ■ 



of acceptable model parameters. 



4.2. 



Q1009+2956 

The absorption-line profile of this system can be repro- 
duced using 3 kinematic components (Fig. f2] and Table [2), 
which are blended together. In addition to N V, there are three 



main constraints: Lya, Cm A977, and S vi A933. We ini- 
tially assume that the three components have the same param- 
eters. We find that using a solar abundance pattern, we can al- 
ways find overlapping regions in parameter space between the 
best fits of Lya and C III A977 if log(n H /cm" 3 ) < 13. These 
solutions are at 2.3 < log (Z/Z ) < 2.5 and 0.2 < logf/ < 0.4, 
(Figs. [7c] andl7dl). However, we cannot find any overlapping 
parameter combinations between suitable fits for S VI X933 
and C III A977. The S VI A933 solutions also do not intersect 
with acceptable Lya solutions, except when log(«H/ cm 3 ) ~ 
12. For most densities, only if we lower the sulphur abun- 
dance by about 1 dex (relative to the solar abundance pattern) 
can we obtain acceptable model profiles for all three of the 
constraints (Table@{. In the case that the absorption at the po- 
sition of S VI A933 is contaminated by blends with unrelated 
lines, an even larger sulphur abundance reduction would be 
required. In principle, it is possible to produce less S VI A933 
by increasing the nitrogen abundance, since N V is the opti- 
mized transition, however we find that such an adjustment is 
not sufficient to fully resolve the discrepancy. For the special 
case of log(« H /cm -3 ) « 12, the S VI A933 and Lya profiles 
can both be reproduced for log (Z/Z ) ^2.3 and \ogU ^0.7. 
For that case, with a higher ionization parameter, an alterna- 
tive acceptable model is found if we raise the carbon abun- 
dance by 0.5 dex, leaving the sulfur abundance unchanged 
from the solar pattern. These results are summarized in Ta- 
bleS 

We next tune separately the metallicity and ionization pa- 
rameter of each component, assuming that all three compo- 
nents have the same volume density, but allowing this com- 
mon value of the density to vary. We select the best model 
parameters at each volume density value using % 2 technique 
and visual inspection (Table 0. An example of acceptable 
models is presented in Fig. [2] To obtain this improved fit we 
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FIG. 7a. — Plots of observed and synthesized profiles of all kinematic 
components of Si IV A 1394 for the absorption system in the spectrum of 
HE0130— 4021. We set log(«H/cm~ 3 ) = 8 and illustrate how the absorption 
profiles change as (log(Z/Z0),logt/) vary about (1.7,-1.5). 
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FIG.7B. — Plots of observed and synthesized profiles of all the kine- 
matic components of Lya for the absorption system in the spectrum of 
HE0130— 4021. We set log(«H/cm~ 3 ) = 8 and illustrate how the absorption 
profiles change as (log (Z/Z : :-,), log £/) vary about (1.7,-1.5). 



vary the metallicities and ionization parameters of different 
components over a range about 0.2 dex. Through this exer- 
cise, we require that S VI A933 is not overproduced. The left 
wing of Lya cannot be well modelled and probably requires 
another component which may not be detected in metal-line 
absorption or may be related to a region to the blue of the 
N V A 1239 profile. We had neglected that in our fit, since we 
could not consider whether it was consistent with N V A 1243. 
Nevertheless, the modelling results for this system indicate 
metallicities considerably higher than the solar values. 



FIG. 7c. — Plots of observed and synthesized profiles of all the kine- 
matic components of Lya for the absorption system in the spectrum of 
Q1009+2956. The different panels show how the models change as 
(log (Z/Z ), log U) vary about (2.5,0.2). In each panel, the model color 
changes gradually from red to blue for logfnn/crn 3 ) = 2, 4, 6, 8, 10, and 
12. 
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FIG. 7d. — Plots of observed and synthesized profiles of all the kine- 
matic components of C III A977 for the absorption system in the spectrum 
of Q1009+2956. The different panels show how the models change as 
(log (Z/Z; ), log U) vary about (2.5,0.2). In each panel, the profiles are color- 
coded according to density in the same way as Fig.l7cl 

4.3. HS1 700+6416 

This absorption system can be described by 3 kinematic 
components (Fig. [3] Table [2]). Because the third compo- 
nent, at v ~ 315kms~', is well separated from the cen- 
tral two components, we determine its parameters indepen- 
dently. Assuming the central two clouds to have the same 
parameters, the preliminary best fits, consistent with both 
Lya and Civ A 1548 as well as with undetected transitions, 
are tabulated in Table [4] In Figures |7e] and [7f] we plot the 
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TABLE 6 

Best model parameters for the HS 1 700+64 1 6 

ABSORPTION SYSTEMS 



Kinem. 


log 




log 


Abundance 


r a 


Comp. 


(z/z ) 


logU 


(« H /cm~ 3 ) 


Pattern 


(pc) 


1 


1.3 


0.5 


2-8 


solar 


9-9000 


2 


1.6 


0.6 




solar 




1 


1.2 


0.5 


10 


solar 


~1 


2 


1.5 


0.6 




solar 




3 


1.7 


1.0 


2-8 


solar 


5-5000 



a The distance of the absorbing gas from the ionizing source for this combination 
of U and n\\ . 

profiles of the two central components of Lya and C IV in 
the log (Z/Z©)-logt/ parameter slice around this preliminary 
best fit for log(n H /cm- 3 ) = 2, 4, 6, 8, 10, and 12. For 2 < 
log(«H/cm~ 3 ) < 8, the line strengths are almost independent 
of density, but for higher densities we see weaker Lya and 
slightly stronger Civ absorption. When log(«H/cm~ 3 ) > 11, 
the acceptable fits to Civ A 1548 occur at log J/ ~ 0.6, but 
Lya is underproduced by the models so that it is impossible 
to obtain an acceptable fit at these very high densities. Ob- 
servationally, variability on a timescale of 6.5 months in the 
quasar's rest-frame implies log(«H/cm -3 ) > 3 (based on con- 
siderations of the recombination time; see 

To refine our preliminary models, we next allow the cen- 
tral two components to have different parameters. Thus, we 
adjust log (Z/Zq) and \ogU of each cloud separately, within 
±0.2 dex of the preliminary best fit. A complete tabulation of 
the best model parameters at each value of the volume density 
can be found in Table[6] Take log(«H/cm -3 ) = 8 for exam- 
ple, a better fit is found at log(Z/Z©) = 1.3, \ogU = 0.5 for 
the first component and log (Z/Z©) = 1.6, \ogU = 0.6 for the 
second component. Even for this example, shown in Fig. [5J 
the Lya components at v ~ km s _1 are slightly underpro- 
duced. This may imply that an additional component would 
be needed, or that the Lya is blended with a line from a sys- 
tem at another redshift. The latter is not unlikely given the 
unidentified component present at v ~ 70 km s 

The best fit for the third kinematic component at v ~ 
315 km s -1 is achieved at log (Z/Z©) = 1.7, logt/ = 1 at 
log(«H/cm~ 3 ) = 2-8. We cannot obtain any acceptable fits 
for log(n H / cm ~ 3 ) > 8 because model Lya profiles become 
broader than the observed profiles for cases that match the 
observed C IV (Figures. |7g"|andl7hii. Combining all three com- 
ponents, Fig.|3]shows one of the acceptable fits we obtain for 
this system. 



4.4. The O VI Doublet 

The energy required to produce Ovi is much higher 
(8.371 Ry) than Civ (3.520 Ry) and Nv (5.694 Ry). Thus, 
the Ovi doublet probes highly-ionized gas and could, in 
fact, arise in a more tenuous gas that is not co-spatial with 
the gas responsible for the lower-ionization absorption lines 
and also has a different coverage fraction of the background 
source(s). With the above in mind, we note that in the case 
of HS 1700+6416 (Fig. |3), the best model cannot fully re- 
produce the O VI absorption in both the central and redward 
components. This could signal that the O VI absorber has a 
larger coverage fraction. Another possibility is that the O VI 
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FIG. 7e. — Plots of observed and synthesized profiles of the two central 
kinematic components of Lya for the absorption system in the spectrum 
of HS 1700+6416. The different panels show how the models change as 
(log (Z/Zp),\ogU) vary about (1.4,0.5). In each panel, the profiles ar e color- 
coded according to density in the same way as Fig. [7c] Figures !7aH7dl |7g| 
and!7hlare available in the online version of the Journal. 
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FIG. 7f. — Plots of observed and synthesized profiles of the two central 
kinematic components of C IV A 1548 for the absorption system in the spec- 
tram of HS 1700+6416. The different panels show how the models change 
as (log(Z/Z.0),log{7) vary about (1.4,0.5). In each panel, the profiles are 
color-coded according to density in the same way as Fig.lTcl 



is contaminated by lines from the Lyman forest, which is sup- 
ported by the large widths of the "black" troughs around the 
Ovi 11038 line of HS 1700+6416 and Q1009+2956 ( Fig- 
ures. |2]and[3]). Almost all spectral regions around O VI for the 
three quasars have some indication of blending, which makes 
the comparison of the synthesized O VI doublet profile to the 
data ambiguous. Therefore, we can only use the observed 
O VI absorption troughs to set an upper limit on the absorp- 
tion produced by our models. 
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FIG. 7g. — Plots of observed and synthesized profiles of the kinematic com- 
ponent at 315 km s of C IV A1548 for the absorption system in the spec- 
tram of HS 1700+64 16. The different panels show how the models change as 
(log (Z/Za),logU) vary about (1.7,1.0). In each panel, the profiles are color- 
coded according to density in the same way as Fig. [7c] The vertical axis scale 
is expanded in order to display this weak feature. 
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FIG. 7h. — Plots of observed synthesized profiles of the kinematic com- 
ponent at 315 km s of Lya for the absorption system in the spectrum of 
HS1700+6416. The velocity zero point is defined as in Fig.[T] In each panel, 
the profiles are color-coded according to density in the same way as Fig. [7c] 
The vertical axis scale is expanded in order to display this weak feature. 

5. PARTIAL COVERAGE OF CONTINUUM AND BROAD-EMISSION 
LINE REGIONS 

In principle, the intrinsic absorbers could cover differ- 
ent fractions of the continuum source and the BELR (e.g., 
Gang uly et alj fl999). Thus an absorption line superposed on 
a broad emission line represents the absorption of different 
fractions of emission-line and continuum photons. If this is 
true, the coverage fraction one obtains from fitting a doublet 
is an effective coverage fraction, which represents the frac- 



tion of the photons from all the background sources that pass 
through the absorber. Thus, for an absorption line that is su- 
perposed on a broad emission line and the continuum, a more 
general expression for the normalized residua l flux at a certain 
velocity is given by (see lGanguly et al.l fl999) 



R(V) 



1 



^ [C c (v) + W(v)C e (v)] 



1+W(v) 



(2) 



(see also th e discussion and more general formalism in 
iGabel. et al.M2.Q05b . In the above expression, C c and C e rep- 
resent the coverage fractions of the continuum source and of 
the BELR and W = fe/fc is the ratio of the flux contribution 
of the emission line (without continuum) and the continuum 
only at the position of the absorption line. In principle C and 
W are functions of velocity across the absorption-line profile. 
However, since the absorption lines of interest here are fairly 
narrow, W does not change appreciably over the profile of the 
absorption line and we can regard it as independent of veloc- 
ity. By rearranging equation (0 one can show that the effec- 
tive coverage fraction, Cf, can be expressed as the weighted 
average of C e and C c as 



C f = 



C c + WC e 
l+W 



(3) 



In practice, we can determine Cf by fitting the profiles of the 
UV resonance doublets and W by measuring the strength of 
the broad emission line relative to the continuum at the loca- 
tion of the absorption doublet. This leads to a relation between 
C e and C c but it does not allow us to determine C e and C c sep- 
arately. An example of a graphical depiction o f this relation 
can be found in Fig. 10 of Gangul y et alJ (119991) . 

The N V absorption doublet presents a further complica- 
tion because it is likely to be superposed on the blend of the 
Lya and N V emission lines. The regions emitting Lya and 
N V are likely to have a different spatial extent and their cov- 
erage fractions should also be different. This is suggested 
by the r esults of reverberation mapping studies (see, for 
example iPeterson & Wandelll2000t lOnken & Peterson! I2001 
iKollatschnvl 12003b . which indicate an ionization stratifica- 
tion in the BELR. Under these circumstances, in the nu- 
merator of equations and (0, we must replace W C e by 
WiyaC^ya + WnvCnv and in the denominator, we must re- 
place W by Wiya + Wnv- In this new notation C\^ a and Cnv 
are the coverage fractions of the Lya and N V emitting regions 
and W^ya and Wnv are defined in a manner analogous to the 
definition of W. 

Putting aside the above complication, we can invert equa- 
tion (0 to obtain an expression for C c , in terms of Cf, W, and 
C e . However, since C e is unknown in practice, this relation 
can only lead to upper and lower bounds on C c , corresponding 
to C e = and 1, respectively. Even if the absorption doublet 
is superposed on a blend of emission lines, the expression for 
Cf reduces to equation (0 if the BELR is fully covered or not 
covered at all. Taking the error bars on the measured quanti- 
ties into account we obtain the following expressions for the 
limits on C c : 



C™ x =min {1, [1 + (W + 5w)] (Q + 8 Cf ) } and (4) 
C c min - max {0, [1 + (W - 5w)} (Q - 8 Cf ) -(W + 5w)} , 
where 5x denotes the error bar on X. 
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TABLE 7 

Constraints on Continuum and BELR Coverage 
Fractions 



Quasar 



Ion 



HE0130-4021 Nv 



LI ±0.4 



Lya 

(C f = 



0.5 ±0.3 



unconstr. 
<0.5 

unconstr. 
<0.4 
<0.3 
<0.4 

unconstr. 
<0.6 
>0.3 

unconstr. 



0.3) 1 



Q1009+2956 Nv 0.66 ±0.12 unconstr. 

>0.2 
<0.5 

Lya 1.2±0.2 < 0.6 unconstr. 

(Q = 0.5) b 

HS 1700+6416 Nv 0.96 ±0.06 0.2-0.8 

< 0.7 
<0.8 

CIV 0.41 ±0.02 0.2-0.8 > 0.4 

Lya 0.6 ±0.2 < 0.7 > 0.2 

(C f = 0.4) b 



a The values of W were obtained from the following spectra: for HE0130— 4021 
we used the spectrum published by Osmer & Smith i 1976). for Q1009+2956 we 
used the spectrum published by Buries & Tvtlei i 1998), and for HS1700+6416 
we used the SDSS DR5 spectrum. 

b In the case of Lya we cannot determine Cf from the data, therefore we assign a 
value equal to the average of all the N V kinematic components. 

In an analogous manner, we can use equation ([3]) to derive 
limits on C e , using the fact that C c must lie between and 1. 
The corresponding expressions, including uncertainties on the 
measured quantities, are: 

, [l±(W + g ty )](C f + g Cf ) | 
15 W^w) J ^ 

\ {W + 5w)-\ J 

We have used flux-calibrated, low-resolution spectra of the 
three quasars to measure the values of W for the Lya, and N V 
absorption lines (as well as C IV absorption line in the case of 
HS 1700+6416). Using these values of W and the values of 
Cf from Table |2] we have estimated the limits on the coverage 
fractions of the continuum source and the BELR, C c and C e , 
which we list in Table[7j Limits derived from the Lya absorp- 
tion line involve the assumption that the effective coverage 
fraction of the Lya absorber is the average of the values mea- 
sured for the N V absorber, which need not be correct. Never- 
theless, these limits appear to be consistent with those derived 
from the N V and C IV doublets. The limits on the continuum 
coverage fraction derived from different N V kinematic com- 
ponents in the same quasar appear to be consistent with each 
other with one exception: component 9 of HE0130— 4021, for 
which the effective coverage fraction is measured to be Cf = 1 . 
The only reliable limit on the BELR coverage fraction comes 
from the C IV doublet of HS 1 700+64 16. 

A closely related issue is whether the superposition of ab- 
sorption lines on strong emission lines dilutes their strengths 



by different amounts and leads to large errors in our estimated 
abundances. The relative strengths of the Lya and Nv ab- 
sorption lines are particularly sensitive to this effect, espe- 
cially when the blueshift of the absorption lines is small. In 
a specific scenario where the absorber covers only the con- 
tinuum source, the Lya absorption line is superposed on the 
peak of the Lya emission line and can be diluted much more 
than the N V line which is superposed on the red wing of the 
Lya emission line (the velocity difference between the Lya 
and N V absorption lines is approximately 5925 km s~'). In 
such a case, the procedure we follow here to determine col- 
umn densities, will lead us to underestimate the column den- 
sity of Lya relative to that of N V and our models would yield 
a nitrogen abundance that is higher than the true one. 

Examining the values of W listed in Table [7j we see that if 
such a scenario is true, the N V absorption line is actually di- 
luted more than the Lya absorption line in HE0130— 4021 
and HS 1700+6416. In other words, in these two quasars 
we would have underestimated the nitrogen abundance rather 
than overestimated it. To verify our expectation quantitatively, 
we explored a new set of models for HS 1700+64 16 in which 
we assumed that the absorber covers the continuum source but 
not the BELR. Thus, we started by adopting the Cf value of 
Cloud 2 because it is the dominant component and set C e = 
and C c =0.64 for all transitions. As a result, the effective 
coverage fraction for each transition depends on its W value, 
e.g., C f = 0.33 for Nv, C f = 0.45 for C iv, and C f = 0.34 for 
Lya. By setting log(«H/cm~ 3 ) = 8, we find an acceptable fit 
at (log (Z/Z ), log U) = (1.6,0.8). Compared to our best fit 
with C c = C e , this metallicity is 0.2 dex higher and the ioniza- 
tion parameter is 0.4 dex higher. 

In Q1009+2956 we have the opposite situation, in which 
WLya is larger than Wnv, i.e., Lya is diluted more by emis- 
sion line than N V, suggesting that metallicity may have been 
overestimated. To assess whether this is indeed the case, 
we explored a new set of models in which the C e has the 
minimum value possible and C c has the maximum value 
possible. According to Eq. and ©, these values are 
C c = 1 and C e = 0.245. Thus, the effective covering fac- 
tors are C f = 0.593 for Lya, and C f = 0.895 for C III. By 
setting log(nn/ cm 3 ) = 8, we find an acceptable model at 
(log (Z/Z ). log U) — (2.5,0.4), very similar to our previous 
best-fitting model which had C c = C e . 

In the calculation above, we assumed the same average 
value of W of two members of the N V doublets. However, the 
value of W could differ between the blue and red members of 
the N V doublet, which could affect not only the metallicity 
estimate but also values of the coverage fractions. In this par- 
ticular case, however, the value of W is the same for the two 
members of the double, within uncertainties. 

6. DISCUSSION 

6.1. High Metallicities 

One of the main results of this work is that the metallicities 
in intrinsic NAL systems are quite high, log(Z/Z Q ) > 1. Sev- 
eral possible scenarios have been proposed to produce super- 
solar metallicities in the vicinity of the cent ral massive black 
holes of quasars (Hamann & Ferlandl [l999l) . The most natu- 
ral scenario involves normal g alactic chemical evolution (see 
lHamann & Ferlandll992l[l993h . By considering spectral syn- 
thesis and chemical enrichment models of the N v/C IV and 
N v/He II emission line ratios, it was suggested that nitrogen 
is over-abundant by factors of 2-9 in high redshift (z > 2) 
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quasars (Hamann & Ferland 1999). Highly evolved gas with 
log(Z/Z ) > 0.5-1.1 is needed to produce such a large over- 
abundance of nitrogen. This evolutionary model requires a 
large number of high mass stars and implies a fiat initial mass 
function (IMF) and rapid star formation. The power-law in- 
dex of the IMF (O °c M T ) can be in the range y = 0.9 - 1 .4 
and star formation must occur in less than 0.5 Gyr for z > 4 
objects (Haman n & Ferlaridll 19921) . This vigorous star forma- 
tion is consistent with that found in the environment of ac- 
tive galactic nuclei (AGNs). The required time scales, metal- 
licities, and IMFs are similar to those of e lliptical galaxies 
and bulges of disk galaxies (Romano et al. 2002). Quasars 
with more massive host galaxies tend to be m ore metal rich, 
based on observations of m assive ellipticals dClemens et al.l 
l2008t iDe Lucia et al.ll2006l) . Therefore, the high abundances 
derived from our absorbing systems are not surprising, as- 
suming this scenario. The quasars could undergo extensive 
chemical evolution before they are observable. Other possi- 
ble processes within the accretion disk that can lead to su- 
persol ar metalicities include metal enrichment by sta r trap- 
ping dArtvmowicz. Lin ~& Wampler 1992; Shields 1996j), star 
formation (Collin & Zahnj |1999tlGoodman & Tanll20o4 . and 
non-stellar nucleosynthesis (Kundt 1996). However, it is un- 
likely that these mechanisms can produce metallicities as high 
aslog(Z/Z )~l. 

Early measure ments of extreme ly high quasar metallicities 
based on BALs dTurnshe k 1986) could have been overesti- 
mated because partial coverage of the background emission 
was not taken in to account, leading to un derestimation of 
column densities (lHamann & Ferlandl[l99 9). However, since 
our analysis of NALs takes partial coverage into account, our 
measurements do not suffer from this bias. Actually, super- 
solar metallicities have also been deriv ed for associated ab- 
sorption systems i n low redshift AGNs dPapovich et ai]l2000t 
Fields et al.ll2005HGabel. Arav & Kiml2006l) . However, those 
metallicities are not higher than log(Z/Z0) ~ 1, i.e., some- 
what lower than our results for quasar intrinsic NAL sys- 
tems. This difference could be due to the higher redshifts 
of our NAL systems and/or due to their presence in more 
luminous AGNs. Rapid and vigorous stellar evolution is 
likely to accompany quasar activity, as suggested by a cor- 
relation between the Eddington ratio and the BELR metal 
abund ance spanning 3 ord ers of magnitude in both quantities 
(e.g., [Shemmer et al. 2004). The resulting high-metallicity 
gas could then be diluted with a low-metallicity interstellar 
medium as it moves outward in the host galaxy. Eventually, it 
can escape the host-galaxy potential, enriching the interstellar 
and intergalactic medium. 

It is also possible that we are looking only at the tip of the 
iceberg, i.e., at the most highly enriched gas. A diversity of 
ab sorber properties , including metallicity, has been reported 
by Petitiean (1996). Not only can absorbers have different 
coverage fractions for the continuum and BELR but also dif- 
ferent ions/species can have different coverage fractions due 
to the multi-phase structure of the gas. Different coverage 
fractions for different ions/species are also inferred from the 
general properties of the strong C IV family of intrinsic NALs 
(Misa wa et al.l 12007). which usually have a saturated Lya 
line, but no detected N V, implying a lower ionization state 
for the gas. In the strong-Nv NAL systems studied in this 
paper, the Ovi doublet may sometimes be under-produced 
by the models, which may be the result of a greater coverage 
fraction of the O VI absorber. Thus the high metallicities that 
we derived for the NAL gas in which N V is detected need not 



apply to all of the gas surrounding the central engine. 

6.2. Properties of the Absorbing Gas 

The ionization parameters that we obtained are similar to 
those for BALs. A typical value of a BAL ionization pa- 
rameter is logf / ~ ( with a range of —1.0 to +0-6; e.g. , 
[ Crenshaw etal] 120031: IWampler, Chugai & Petitieanl 1 19951: 
iTelfer et al.l 119981) . Both an intrinsic NAL and a 
BAL are observed in the spectru m of Q0059— 2735 by 
Wamrjler, Chugai & Petiti eanl (Il995[). who s how th at both can 
be fitted with logf/ = -0.7 iHamamil d 19981) finds that 
logf/ > -0.6 is required to model the BAL of PG 1254+047, 
while lTelfer et al.l 0998), using a two-slab geometry to model 
the BAL system of QSO SBS 1542+541, find -1.1 < logf/ < 
0.6 for the lower-ionization zone and logf/ > 0.3 for the 
higher-ionization zone. 

If we assume that the outflow includes thin filaments (jus- 
tified below) that are responsible for the UV absorption lines 
that we observe, we can estimate the outflow rate in these fil- 
aments as 

where p is the typical mass density of absorbing filaments, 
r is their typical distance from the center, 8r is their typical 
thickness in the radial direction and SQ. is the solid angle they 
subtend to the center of the flow. We can re-cast the above 
expression by using the definition of the ionization parameter, 
U = Q/4-nr 2 nnc, to replace 4-nr 2 p (Q is the rate of emission 
of ionizing photons at E > 1 Ry, and «h is the hydrogen num- 
ber density). We also note that the thickness of a filament 
is related to the hydrogen density, and radial column density, 
Nn, via A^h = n dr. Thus we can write 

M = 5.2xlO- 5 M 8 U 2 /_ 1 -^f- Mgyr- 1 , (8) 

\n 6 U J 

where A^h = 10 18 Ni 8 cm" 2 , v = I0 2 v 2 kms" 1 , (8Q./4n) = 

10 _1 /-i. Q = 10 58 658 s _1 , andn H = 10 6 n 6 cm" 3 . We also 
note that from the definition of the ionization parameter, we 
can obtain an estimate of the distance of the absorbing fila- 
ment from the ionizing source. 

To justify the assumption that the absorbers are thin, we 
note that using the definition of U and A^h = n8r, we 
obtain 8r/r = N^AkcU /n\\Q. The values of Nu (neu- 
tral plus ionized) inferred from our photoionization mod- 
els are in the following range : log(AVcm- 2 ) ~ 16-19 for 
HE01 30-4021, log(Ar H /cm- 2 ) ~ 17-18 for Q1009+2956, 
and log(N H /crrr 2 ) - 19 for HS 1700+6416, the values of U 
are given in Tables[4]-[6]and the values of Q are given in Table 
[TJ Thus, any combination of values of the above quantities 
yields 8r/r < 0.01 which indicates that the absorbers are ge- 
ometrically thin along the line of sight. 

We can now use equation [8] to obtain constraints on the 
mass outflow rate. The necessary values of A^h, V, U and 
Q are given in Tables [4]-[6] We also use a v alue of /_i ^5 
based on the results of Misawa et al. (2007). Unfortunately, 
we have very poor constraints on the hydrogen volume den- 
sity, which lead to a very wide range of values for the mass 
outflow rate, i.e. M ~ few M e yr~ ! for the lowest densities 
and M < 10~ 4 M yr~' for n H > 10 10 cm" 3 . The kinetic 
power of these filaments is negligible compared to the elec- 
tromagnetic luminosity (£k = jMu 2 /Lboi < 10~ 5 ), however 
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they may be embedded in a hotter, more massive outflow or 
they may be accompanied by a more massive outflow along a 
different line of sight. 

Finally, we note that the poor constraints on the density also 
lead to poor constraints on the distance of the filament from 
the ionizing source, r. In Tables [5]-|6] we list the constraints 
on the radial distance of the filaments as a function of density 
for Q 1009+2956 and HS 1700+64 16 (the results for HE0130- 
4021 are very similar). For«H > 10 9 cm~ 3 , the absorbers turn 
out to be very close to the ionizing source, at distances com- 
parable to the size of the BELR, while at the lowest densities, 
the filaments are located in the outskirts of the host galaxy. 

7. SUMMARY AND CONCLUSIONS 

In this paper, we modelled the intrinsic NAL systems in 
the spectra of the quasars HE0130-4021, Q1009+2956, and 
HS 1700+6416, We identify these systems as intrinsic be- 
cause they exhibit partial coverage in the N V AA1239, 1242 
doublet. Using the photoionization code Cloudy, we simulate 
and reconstruct each absorption feature observed in a system 
by adjusting the metallicity (relative to the Sun), Z/Z , ion- 
ization parameter, U and hydrogen volume density, «h- Our 
main conclusions are: 

1. All three intrinsic systems have metallicities of 
log(Z/Z Q ) > 1-2, regardless of constraints on U or «h- 
We find the high metallicities that we infer cannot al- 
ternatively be explained by different coverage fractions 
of the continuum source and BELR. The origin of these 
supersolar metallicities is unknown, but we speculate 
that they may result from intensive star formation be- 
fore the quasar becomes observable. They could also 
be related to an inhomogeneous metallicity distribution 
in the vicinity of the quasar central engine. 

2. All three intrinsic systems require high ionization pa- 
rameters, log U ~ 0, which are similar to those derived 
for BALs. 

3. We can constrain the coverage fractions of the contin- 
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uum source and the BELR separately, and find that the 
continuum source is not fully covered by the absorbers. 

4. Although we cannot constrain the hydrogen volume 
density very well, we do find that the widths of the syn- 
thesized line profiles become larger than the observed 
ones at very high densities. This is a result of a sharp 
increase in temperature at high densities, produced by 
free-free absorption of infrared photons by free elec- 
trons. This leads us to conclude that the hydrogen vol- 
ume density is in the range 2 < log(«H/cm -3 ) < 9. 

5. As a result of the large range of possible densities we 
cannot obtain robust estimates of the mass outflow rate 
in the filaments responsible for the N V absorption. We 
can conclude, nevertheless, thatM < few M yr , and 
the corresponding power is a negligible fraction of the 
electromagnetic luminosity. 

6. The strength of the O VI AA1032, 1038 doublet is un- 
derproduced by the models that reproduce the other ab- 
sorption lines. This suggests that the O VI lines are ei- 
ther contaminated by lines in the Lyman forest of that 
they arise in a separate gas phase with a larger coverage 
fraction. 

The strong-Nv absorbers studied in this paper make up 
one of two, a pparently distinct, fa milies of intrinsic absorber 
identified by (Misawa et al. 2007). Our next goal is to study 
the physical properties of the other family of absorbers, the 
strong-C IV family and probe the origin of the two families. 
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